Clocking hadronization in relativistic heavy ion collisions with balance functions 
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A novel state of matter has been hypothesized to exist during the early stage of relativistic heavy 
-- - ' ion collisions, with normal hadrons not appearing until several fm/c after the start of the reaction. 

' To test this hypothesis, correlations between charges and their associated anticharges are evaluated 

, with the use of balance functions. It is shown that late-stage hadronization is characterized by 

' tightly correlated charge-anticharge pairs when measured as a function of relative rapidity. 

Relativistic heavy ion collisions produce mesoscopic regions of enormous energy density, perhaps surpassing 3 
^ ■ GeV/fm^ in Pb collisions at the CERN SPS with even higher energy densities expected at RHIC. At such energies 
I hadronic degrees of freedom should be replaced by quark-gluon degrees of freedom. Several experimental measurements 
• have been proposed as signals to the quark-gluon plasma Among these signals is an expected enhancement in 
strange-quark production which should take place 5-10 fni/c into the collision when the local temperature has dropped 
to near 160 MeV, but the system is still far from freeze-out. Strangeness enhancement has indeed been observed 
OA in heavy ion collisions but alternative hadronic explanations have also been put forward assuming early-stage 
^ , hadronization with medium modifications, referred to as color ropes |^,^ or baryon junctions |^]. In this paper the 
N ' use of balance functions is proposed as a means to determine whether quark production occurred at early times, t < 1 
i fm/c, or according to a late-stage hadronization scenario, see e.g. [||,^. 

' Late-stage production of quarks could be attributed to three mechanisms: formation of hadrons from gluons, conver- 
' sion of the non-perturbative vacuum energy into particles, or hadronization of a quark gas at constant temperature. 
■ Hadronization of a quark gas should approximately conserve the net number of particles due to the constraint of 
' entropy conservation. Since hadrons are formed of two or more quarks, creation of quark-antiquark pairs should 
, accompany hadronization. All three mechanisms for late-stage quark production involve a change in the degrees of 
freedom. Therefore, any signal that pinpoints the time where quarks first appear in a collision would provide valuable 
insight into understanding whether a novel state of matter has been formed and persisted for a substantial time. The 
fact that the hadronic phase has a higher concentration of charges than the QGP phase at the same entropy has been 
^ , discussed in the context of charge fluctuations in ||l^ . 
• • ■ The link between balance functions and the time at which quarks are created has a simple physical explanation. 
Charge-anticharge pairs are created at the same location in space-time, and are correlated in rapidity due to the 
strong collective expansion inherent to a relativistic heavy ion collision. Pairs created earlier can separate further in 
rapidity due to the higher initial temperature and due to the diffusive interactions with other particles. The balance 
5^ 1 function, which describes the momentum of the accompanying antiparticle, quantifies this correlation. 

The balance functions employed here are similar to observables used to investigate hadronization in jets produced 
in pp or e+e~ collisions |l^Jl^. The balance function describes the conditional probability that a particle in the bin 
Pi will be accompanied by a particle of opposite charge in the bin p2. We define the balance function, 

B{p2\pi) = ^ {pib,p2\a,pi) - p{b,p2\b,pi) + p{a,p2\b,pi) - p{a,p2\a,pi)} , (1) 

where p(b,p2\a,pi) is the conditional probability of observing a particle of type b in bin p2 given the existence of 
a particle of type a in bin pi. The label a might refer to all negative kaons with b referring to all positive kaons, 
or a might refer to all hadrons with a strange quark while b refers to all hadrons with an antistrange quark. The 
conditional probability p{b,p2\a,pi) is generated by first counting the number N{b,p2\a,pi) of pairs that satisfy both 
criteria and dividing by the number N(a,pi) of particles of type a that satisfy the first criteria. 



X N{b,p2\a,pi) 

p{b,p2,a,pi) ^ — (2 

Both sums run over all events, though pairs only involve particles from the same event. 

An example of binning might be that pi refers to a measurement anywhere in the detector, while p2 refers to 
the relative rapidity \yi, — ya\. Then the balance function would be a function of Ay only, and would represent the 
probability that the balancing charges were separated by Ay (in our formalism we include a division by Ay to express 
B{Ay) as a density). 

The balance function is normalized to unity if a/b refer to all particles with a positive/negative globally conserved 
charge. 



B{P2 \Pi)^l {Mb - (Mb - 1) + Ma - {Ma - 1)} = 1, (3) 

P2 

where Ala and Mb are the average multiphcities of the a and b particles. The normahzation derives from the fact that 
for every extra positive charge there exists one extra negative charge. If the acceptance measures only a fraction of 
the charge, e.g. only kaons are measured and the strangeness in hyperons is excluded, the balance function would sum 
to that fraction. Balance functions can exploit any conserved charge: electric charge, strangeness, baryon number or 
charm. The first two terms in Eq. constitute the balance functions defined in several analyses of e~^e~ —> jets. 
By adding the last two terms the normalization properties are retained even for the case where there is a non-zero 
net charge, Ala — AI^ ^ 0. 

If many charges are present in the event, the balance function represents the subtraction of two large numbers. 
However, large multiplicities also imply a large number of pairs from which to calculate the balance function. Since 
the number of uncorrelated pairs rises as the square of the multiplicity Af , the statistical error in calculating the 
numerators of the conditional probabilities, which rises as the square root of the number of pairs, increases linearly 
with M . Since the denominator also rises linearly with Af, the statistical error in the balance function is independent 
of multiplicity and is principally determined by the number of events: 
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Thus, the baryon-antibaryon balance function which might involve a few dozen antibaryons would require the same 
number of events as the electric-charge balance function which might be constructed from a thousand particles. 
Typically, 10^ events are required to determine a balance function with statistical fluctuations at the level of 10^^. 

Balance functions probe the dynamics of charge-anticharge pairs by quantifying the degree to which the charges 
are correlated in momentum space given the constraint of being created at the same space-time point in a system 
exhibiting strong position-momentum correlations such as a relativistic collision where source velocities might span 
several units of rapidity. In a globally equilibrated system with no collective flow, there would exist no correlation 
between the balancing charges, and the numerator in Eq. would factorize. The width of the balance function 
would then correspond to the extent of single-particle emission in momentum space. 

To illustrate the way in which balance functions quantify the charge-anticharge correlations, we consider a Bjorken 
boost-invariant parameterization [|l3| of a source expanding along the z axis with a collective velocity proportional 
to the position, Wcoii = zjt. All intrinsic variables, such as density or temperature, depend only on the proper time 
T = ii? — z^)^!"^ . We first consider only direct production of hadrons, as the possibility of hadrons coalescing from 
quarks is discussed later in the paper. Particles and antiparticles of mass m are generated in pairs at the same point 
in space-time following a local thermal distribution, and the relative rapidities are used to generate balance functions. 
The characteristic width of the balance function is determined by the ratio of the temperature to the mass. Non- 
relativistically, Oy = {2T /mY/'^ ^ and heavier particles are characterized by narrower balance functions. For particles 
with masses much less than the temperature, the balance functions become independent of the temperature. 

Figure ^ displays balance functions assuming a Bjorken parameterization of an expanding pion gas and an expanding 
proton gas, for two temperatures, 225 MeV and 165 MeV. Clearly, the balance functions of the more massive particles 
are sensitive to the temperature. This suggests that the strangeness and baryon balance functions should provide 
more insight than the electric-charge balance function which would be largely dominated by pions. 

Balance functions in heavy ion collisions should be compared to those from pp collisions at the same -y/s where 
hadronization is nearly instantaneous. Charged-pion balances measured in e^e~ collisions as a function of the rapidity 
defined along the jet axis have been reasonably explained by the string hadronization dynamics of the Lund model 
iQ, e.g. as implemented in PYTHIA |l^. Thermally generated balance functions are compared to predictions of 
PYTHIA for pp collisions at = 200 GeV in Fig. |l|. The PYTHIA balance functions tend to be broader than 
those that are thermally generated, especially for the more massive protons and kaons. Assuming that experimental 
balance functions in pp collisions would be well described by similar string dynamics. Fig. |^ suggests that narrower 
balance functions might indeed point to thermal production at a lower temperature and thus at later times in the 
evolution of the heavy ion reaction. 

Rescattering and annihilation should also affect balance functions. Rescattering may be qualitatively understood 
by considering the diffusion equation in Bjorken coordinates r and rj = tanh~^(z/t), where 77 plays the role of the 
position in the z direction and also equals the collective rapidity of the local matter. Rather than considering the 
diffusion constant D = vt/{na) as a constant, it is more physical to incorporate the fact that the density n falls 
inversely with r and to consider (3 = vt/ {nra) as a constant where Vt is the thermal velocity and cr is a characteristic 
cross section. The diffusion equation then becomes 



Here, / is the probability of observing a particle at position rj at time r. With the initial condition of ry = at tq, 
the solution to the diffusion equation is a Gaussian with variance cr^ = 2/31n(r/To). This illustrates that collisions 
broaden the balance function by diffusing the charge in the effective spatial coordinate 77. However, in the limit of 
zero mean free path, the diffusion constant tends to zero and the particles do not then diffuse. 

The overall width of the balance function in relative rapidity is a combination of the thermal rapidity spread crthorm 
and the effect of diffusion in 77 of both particles: 

='^thcrm + 4/31n(T/ro). (6) 

Due to cooling, the width CTthcrm falls with time which provides a competition between diffusion which stretches the 
balance function, and cooling which narrows it. If the production occurs at early times, then ln(r/To) is large and 
the effect of collisions is to significantly broaden the balance function. 

Some hadrons will contain coalesced quarks that were created at early times. The thermal contribution to ay 
described in Eq. should be unaffected by the past history of the constituent quarks. However, the diffusive contri- 
bution might significantly depend on the fact that the charge moved as a free quark rather than as a hadron during 
it's early history. Balance functions constructed from hadrons can thus provide meaningful information regarding the 
creation and mobility of the constituent quarks. 

To quantitatively illustrate the effect of rescattering, we model a pair of particles produced at an initial proper time 
To that collide A^coU times before disassociating at a final time r/ . Each collision is assumed to completely reorient the 
particle with the local collective velocity. The collision times are chosen randomly such that the number of collisions 
as a function of ln(T) is uniform. The temperature is chosen to vary linearly with the proper time, cooling from 225 
MeV at T = f fm/c to f20 MeV at r = 15 fm/c. Figure || shows the Kj^K^ balance function with iVcoU = and 
-^coU = 10 assuming kaons are created at r = 1 fm/c and cease to collide at r/ = 15fm/c. In this case collisions 
clearly broaden the balance function. 

Annihilations should also broaden the balance function. Annihilation forms new correlated pairs with the surviving 
partners of the annihilated particles, which tend to be less correlated than the original pairs. Annihilation combined 
with an equal amount of creation does not affect the balance function since the relative rapidities of formed and 
annihilated pairs should be identical. Figure |2| illustrates the effects of annihilation by considering the same case 
described above, but with the additional assumption that half the particles disappear due to annihilation. In hadronic 
models of heavy ion collisions, the number of both antibaryons and strange particles tend to decrease with time due 
to cooling, which should result in broadened balance functions. 

Figure ^ displays the effect of collisions on balance functions for pions, kaons and protons, by considering the 
mean relative rapidity as a function of the number of collisions. For production at early times when the collective 
velocity gradient is high [dvcoii/ dz — l/r), collisions broaden the balance function. However, for very large numbers 
of collisions, the charge does not diffuse and the balance functions are narrowed due to the cooling. One would expect 
particles to undergo 10-20 collisions if created at r = 1 fm/c, although the effective number of completely randomizing 
collisions might be closer to a half dozen. If created at r = 9 fm/c when the temperature is 165 MeV, the effective 
number of completely randomizing collisions might be two or three. Figure ^ suggests that the signal for late-stage 
quark production is significantly magnified by rescattering. Due to collisions, even charged-pion balance functions 
become strongly sensitive to the creation time. 

The simple calculations presented here sidestep two issues: correlations from decays such as </> ^ K^K^, and 
experimental acceptance problems. Both problems can be addressed by modeling constrained by the multitude of 
other observables measured in a heavy ion collision. Although some open questions remain, it seems clear that the 
canonical picture of a heavy-ion reaction,quark-gluon plasma formation followed by late-stage hadronization, should 
have a clear signature in the balance functions. Compared to pp collisions, one expects the peak in the balance 
function in nucleus-nucleus collisions to be narrower near Ay — due to the contribution of late-stage production of 
quark pairs, while the tails of balance function should become broader reflecting the extra diffusion of charge in the 
early stages of the collision. Finally, we remark that we have barely explored the possibilities of balance functions. 
The rich nature of the binnings {p2\pi) should provide a powerful means for resolving many of the issue regarding 
creation and diffusion of quarks and hadrons in relativistic heavy ion collisions. 
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FIG. 1. Balance functions as predicted in a simple Bjorken thermal model are shown for two temperatures, 225 MeV and 
165 MeV. Since heavier particles from cooler systems have smaller thermal velocities, they are more strongly correlated in 
rapidity and result in narrower balance functions. Also shown are balance functions as predicted by PYTHIA where the shape 
of the balance function is largely determined by string phenomenology. 
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FIG. 2. Kaon balance functions arc shown assuming an initial local temperature of 225 MeV and a production time of 1 
fm/c. The balance function is broadened by the inclusion of randomizing collisions and annihilation. 
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FIG. 3. The mean width of the balance function displayed as a function of the number of collisions, both for the case where 
particles are created early (r = 1 fm/c, T = 225 MeV) and late (r = 9 fm/c, T = 165 MeV). 



